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Abstract:  This  final  technical  report  begins  with  the  objective,  followed  by  detailed  reports  on  new 
accomplishments  and  new  findings  obtained  in  the  project.  They  are  presented  in  six  sections:  A. 
Potential  Energy  Surfaces  of  Ion-Molecule  Reactions,  B.  Potential  Energy  Surfaces  of  Ions,  C.  Potential 
Energy  Surfaces  for  Photochemical  Reactions  and  Nonadiabatic  Processes,  D.  Potential  Energy  Surfaces 
of  Open-Shell  Molecules,  E.  Dynamics  and  Kinetics  in  the  Ground  Electronic  State  and  F.  Potential 
Energy  Surfaces  for  the  Dynamics  of  Chemical  Laser  Systems.  .  The  Publication  List  and 
Interactions/Transactions  conclude  the  report. 


1.  Objectives. 

Understanding  the  mechanism,  kinetics  and  dynamics  of  elementary  gas  phase  reactions  is  one  of 
the  major  goals  of  chemistry.  Such  understanding  is  also  essential  to  predicting  and  understanding 
plasma  dynamics  and  optical  radiation  associated  with  the  spacecraft-atmosphere  interactions.  The 
objective  of  the  present  research  under  the  AFOSR  grant  (F49620-98- 1-0063)  is  to  provide,  based  on 
mainly  the  ab  initio  molecular  orbital  and  some  dynamics  calculations,  theoretical  information  concerning 
the  potential  energy  surfaces  (PESs)  that  dictate  the  kinetics  and  dynamics  of  gas  phase  elementary 
reactions.  The  major  theoretical  thrust  of  the  present  approach  is  to  study  more  than  one  potential  energy 
surfaces  at  the  same  time,  paying  special  attention  to  the  seam  of  crossing  of  potential  surfaces  (conical 
intersection  between  the  states  of  the  same  spin,  or  crossing  between  the  states  of  different  spin)  and  to 
the  non-adiabatic  transition  between  them.  Many  of  the  systems  for  which  we  perform  theoretical 
calculations  are  relevant  to  upper  and  lower  atmospheric  chemistry  and  chemical  lasers,  problems  of 
major  interest  of  the  Air  Force.  They  are  chosen  based  on  the  excellent  experimental  studies  by  groups  at 
Air  Force  Research  Laboratory,  in  order  to  provide  them  with  some  new  insight  that  is  not  easily 
available  without  theoretical  studies. 

2.  Accomplishments  and  New  Findings. 

A.  Potential  Energy  Surfaces  of  Ion-Molecule  Reactions 
A-l.  Mechanism  of  the  C2H2+  +  NH3  reaction. 

Q.  Gui  and  K.  Morokuma,  Ab  initio  study  on  the  mechanism  of  C2H2+  +  NH3  reaction.  Efficient  charge 
transfer  and  proton  transfer  processes  competing  with  stable  complex  formation.  J.  Chem.  Phys.,  108, 
4021-4030  (1998). 

High  level  ab  initio  calculations  have  been  performed  to  investigate  the  mechanism  of  the  ion- 
molecule  reaction  NH3+C2H2+.  Three  channels,  covalent  complex  formation  (CC),  proton  transfer  (PT) 
and  charge  transfer  (CT)  have  been  studied.  Among  the  two  pathways  found  for  the  PT  channel,  one 
leads  the  reactants  NH3+C2H2+  to  NH4++C2H(2n)  through  a  moderately  bound  complex  without  any 
barrier,  and  the  other  leads  NH3++C2H2  to  the  H-atom  transferred  products  NH,++C2H(2X*)  with  a 
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modest  barrier.  These  findings  support  the  fast  “stripping”  mechanism  proposed  by  Anderson  et  al.  As  to 
the  CC  channel,  several  isomers  of  C2H5N+  and  the  isomerization  transition  states  have  been  located.  No 
significant  barrier  relative  to  the  reactants  have  been  found  on  either  the  ground  or  the  2A”  excited  state. 
To  rationalize  the  experimental  fact  that  no  CC  channel  products  have  been  observed,  it  is  argued  that  the 
reactants  NH3+C2H2+  correlate  adiabatically  to  excited  states  of  covalent  C2H5N+  species,  whose 
formation  requires  significant  alternation  of  the  C2H2+  geometry  and  electronic  structure.  Therefore,  the 
system  is  most  likely  to  follow  the  PT  or  the  CT  channel  instead  of  visiting  the  CC  channel.  For  the  CT 
channel,  limited  potential  energy  surface  scans  of  the  three  electronic  states  (1,2  2 A'  +  2 A")  indicate  that 

CT  at  different  approach  angles  or  between  electronic  states  of  different  symmetries  (A’—>A\  A”— »A’) 

may  produce  final  products  of  different  characteristics,  and  might  account  for  the  two  pathways  proposed 
by  Anderson  et  al. 

A-2.  Mechanism  of  the  C2H2+  +  CH4  reaction. 

Q.  Cui,  Z.  Liu  and  K.  Morokuma,  Theoretical  study  on  the  mechanism  of  CH4  +  C2H2+  reaction.  Mode- 
enhancement  effect,  J.  Chem.  Phys.,  109,  56-62  (1998). 

High  level  ab  initio  calculations  have  been  performed  to  investigate  the  mechanism  of  the  ion- 
molecule  reaction  of  CH4  +  C2H2+.  Except  for  some  subtle  differences,  the  profile  for  the  H-abstraction 
channel  obtained  here  at  the  G2M//B3PW9 1/6-31  lG(d,p)  level  is  very  similar  to  that  found  in  a  previous 
study  at  the  G2//MP2/6-31G(d)  level.  For  the  complex  formation  channel,  however,  a  different  transition 
state  has  been  located;  the  geometry  and  energetics  of  which  are  more  consistent  with  experimental 
findings.  Calculations  of  a  few  direct  trajectories  have  been  carried  out  to  investigate  the  possible  reason 
for  the  significant  mode  enhancement  observed  experimentally  for  the  H-abstraction  channel.  Although 
none  of  them  is  reactive,  a  trajectory  with  an  asymmetric  C2H  bend  excitation  exhibits  a  clear  signature 
for  being  more  reactive  than  those  without  vibrational  excitation  or  with  a  symmetric  bend  excitation. 

K.  Morokuma,  Q.  Cui,  and  Z.  Liu,  Potential  Energy  Surfaces  and  Their  Dynamic  Implications  —  Non- 
adiabatic  transition  state  theory  for  spin-forbidden  reaction  CH(2P)  +  N2  and  channel  and  mode  selectivity 
in  ion-molecule  reactions:  C2H2+  +  NH3  and  C2H2+  +  CH4,  Faraday  Discussion,  110,  71-89  (1998). 

Accurate  density  functional  and  ab  initio  calculations  have  been  performed  to  study  the  potential 
energy  surfaces  and  their  implications  to  kinetics  and  dynamics  for  the  following  reactions.  I.  Spin- 

forbidden  reaction  CH(2P)  +  N2  — »  HCN  +  N(4S).  Potential  energy  surface  characteristics  are  calculated 

and  used  to  evaluate  the  overall  rate  using  the  non-adiabatic  transition  state  theory.  II.  Gas-phase  ion- 
molecule  reactions.  A.  C2H2+  +  NH3.  Potential  energy  surfaces  have  been  calculated  and  the  mechanism 
of  efficient  charge  transfer  and  proton  transfer  competing  with  stable  complex  formation  is  discussed.  B . 
C2H2+  +  CH4.  Mode-enhancement  effect  has  been  elucidated  in  terms  of  the  new  transition  state  and  by 
direct  trajectory  calculations. 
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A-3.  Mechanism  of  the  C2H2+  +  methanol  reaction. 

S.  Irle  and  K.  Morokuma,  Ab  initio  and  density  functional  study  on  the  mechanism  of  the  C2H2+  + 
methanol  reaction,  J.  Chem.  Phys.,  Ill,  3978-3988  (1999). 

High  level  ab  initio  (G2MS  and  CASSCF)  and  density  functional  (B3LYP)  calculations  were 
carried  out  to  study  the  mechanism  of  the  ion-molecule  reaction  C2H2+  +  CH3OH  for  four  reaction 
channels:  hydride  abstraction  from  methanol  (HA),  proton  transfer  from  acetylene  cation  (PT),  charge 
transfer  (CT)  and  covalent  complex  formation  (CC)  channel.  For  the  CT  channel,  two  pathways  have 
been  found:  a  usual  non-adiabatic  pathway  via  A7A”  seam  of  crossing  and  a  low  energy  adiabatic 
pathway  through  an  initial  intermediate;  the  latter  may  be  the  dominant  process  with  favorable  energies 
and  a  large  impact  parameter.  The  HA  process  involves  a  low-energy  direct  intermediate  and  a  very  low 
barrier  to  form  C2H3  +  CH2OH+  and  is  also  energetically  favorable.  The  PT  processes  require  passage 
over  a  high-energy  transition  state  and  are  not  important.  One  of  the  experimentally  unobserved  CC 
channels,  formation  of  the  COCC  skeleton,  is  energetically  favorable  and  there  is  no  energetic  reason  for 
it  not  to  take  place;  a  “dynamic  bottleneck”  argument  may  have  to  be  invoked  to  explain  the  experiment. 
The  increase  in  reaction  efficiency  with  the  CC  stretch  excitation  may  be  justified  by  considering  the  TSs 
for  two  CT  pathways,  where  the  CC  distance  changed  substantially  from  that  in  the  reactant  C2H2+.  Very 
qualitatively,  the  C2H2+  +  CH3OH  potential  energy  surface  looks  more  like  that  of  the  C2H2+  +  NH3 
system  than  the  C2H2+  +  CH4  system,  because  of  the  differences  in  the  ionization  potentials:  NH3  ~ 

ch3oh  <  c2h2<  ch4. 

A-4.  The  mechanism  of  the  ion-molecule  reaction  of  0+  (4S)  +  C2H2. 

K.  Fukuzawa,  T.  Matsushita,  K.  Morokuma,  D.  J.  Levandier,  Y.  Chiu,  R.  A.  Dressier,  E.  Murad,  A. 
Midey,  S.  Williams,  and  A.  A.  Viggiano,  An  ab  initio  and  experimental  study  of  vibrational  effects  in 
low  energy  0+  +  C2H2  charge-transfer  collisions,  J.  Chem.  Phys.,  115,  3184-3194  (2001). 

Theoretical  and  experimental  studies  are  performed  to  elucidate  the  low  energy  charge-transfer 

dynamics  of  the  reaction,  0+  (4S)  +  C2H2  (X  'Yg+)  -»  O  +  C2H2+.  In  particular,  the  role  of  the  low- 

frequency  acetylene  bending  modes  (612  and  730  cm'1)  in  promoting  charge  transfer  was  examined. 
High  temperature  guided-ion  beam  measurements  are  carried  out  over  the  energy  range  from  near-thermal 
to  3  eV  at  310  and  610  K.  The  charge-transfer  cross  sections  are  found  to  decrease  up  to  0.5  eV,  to  have 
a  constant  value  at  intermediate  energies  between  0.5  and  1.5  eV,  and  then  to  dramatically  increase  above 
a  threshold  of  a  spin-allowed  process  determined  to  be  at  1.7  eV.  A  bending  vibrational  enhancement  of 
~8  is  observed  at  intermediate  energies.  Thermal  energy  rate  coefficients  are  measured  in  a  variable 
temperature-selected  ion  flow  drift  tube  apparatus  from  193  to  500  K.  At  each  temperature,  a  negative 
energy  dependence  is  observed.  In  order  to  elucidate  the  reaction  mechanism  in  detail,  high  level  ab  initio 
calculations  using  Complete  Active  Space  Self-Consistent  Field  (CASSCF)  and  Multi-Reference  Single 
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and  Double  excitation  Configuration  Interaction  (MRSDCI)  methods  have  been  performed.  The  results 
indicate  that  the  charge  transfer  reaction  occurs  at  an  early  stage  via  nonadiabatic  transition  between 
quartet  and  doublet  states.  There  is  a  weak  van  der  Waals  minimum  at  the  entrance  channel  between  0+ 
(4S)  and  C2H2  with  the  relative  energy  of  -1.51  kcal/mol.  The  minimum  of  the  quartet/doublet  crossing 
seam  (Q/D  MSX),  where  the  spin-forbidden  nonadiabatic  transition  is  most  likely  to  take  place,  lies  very 

near  this  minimum  at  Rc0  =4.06A,  Rcc  =  1.20A,  and  ZCCH  =  166.6°  with  a  relative  energy  of  -1.48 
kcal/mol.  After  the  nonadiabatic  transition  the  system  propagates  on  the  doublet  surface  to  reach  the 
exothermic  O('D)  +  C2H2+(x2riu)  products.  No  energy  barrier  exists  on  the  reaction  pathway,  strongly 

suggesting  that  the  reaction  should  occur  at  low  energy  with  a  negative  energy  dependence,  which  is 
consistent  with  the  experiment.  The  Q/D  MSX  has  a  bent  acetylene  moiety,  which  suggests  that  the 
excitation  in  bending  modes  will  enhance  the  reaction,  in  agreement  with  the  experiment. 

A-5.  Potential  energy  profile  and  mechanism  of  ion-molecule  reaction:  CH4  +  02+  (2ng). 

S.  Irle  and  K.  Morokuma,  Ab  initio  investigation  of  the  potential  energy  profiles  of  gas  phase  CH4  +  02+ 
(2ng)  reaction  system,  J.  Chem.  Phys.,  114,  6119-6127  (2001). 

High  level  ab  initio  and  density  functional  calculations  of  the  ground  state  potential  energy 

profiles  were  carried  out  to  study  the  mechanism  of  the  ion-molecule  reaction  CH4  +  02+  (X  2flg)  for  four 

reaction  channels:  insertion  of  02+  into  the  C-H  bond  of  methane  (INS),  hydride  abstraction  from 
methane  (HA),  charge  transfer  (CT),  and  0-0  cleavage  path  (OO)  after  INS  process  to  give  CH2OH+  + 
OH.  Common  to  these  channels  are  initial  encounter  complexes,  and  our  calculations  match  very  closely 
experimental  estimates  for  binding  energies.  The  INS  channel  proceeds  through  CH400+  and  gives  a 
deep  minimum  corresponding  to  exothermic  and  metastable  intermediate  CH3OOH+.  This  species  can 
easily  eliminate  H  to  give  CH2OOH+,  a  product  observed  experimentally.  For  the  slightly  endothermic 
HA  channel  to  give  CH3+  +  OOH,  two  pathways  were  found:  a  direct  pathway  (likely  to  dominate  at 
higher  collision  energy)  from  the  encounter  complex  via  an  HA  transition  state  at  5.8  kcal/mol  above  the 
reactants,  and  an  indirect  pathway  with  a  slightly  smaller  energy  requirement  consisting  of  elimination  of 
OOH  from  the  INS  intermediate  CH3OOH+.  A  transition  state  with  a  high  energy  requirement  of  15 
kcal/mol  was  found  for  0-0  cleavage  from  CH3OOH+,  consistent  with  the  experimental  finding  that  0-0 
cleavage  occurs  at  high  energies.  It  was  also  found  that  the  seam  of  crossing  between  two  potential 
surfaces  is  facilitated  and  therefore  the  CT  channel  is  promoted  by  the  0-0  stretching  and  the  methane 
deformation  vibrations,  again  consistent  with  the  experimental  results. 

B.  Potential  Energy  Surfaces  of  Ions 
B-l.  Structure  and  spectroscopy  of  03‘  ion. 

Q.  Cui  and  K.  Morokuma,  Ab  initio  studies  on  the  electronic  excited  states  and  photodissociation  of  03 
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anion,  J.  Chem.  Phys.,  108,  7684-7694  (1998). 

Extensive  ab  initio  calculations  have  been  performed  for  the  low  lying  electronic  states  of  Of  to 
elucidate  the  mechanism  of  photodissociation  processes.  The  identity  of  the  mysterious  state  implied  by 
the  recent  experiment  of  Continetti  et  al.  has  been  discussed  based  on  the  current  calculations. 

Calculations  reveals  that  2B2  is  a  truly  bound  state  favoring  a  strongly  bent  geometry  with  ZO-O- 

0-90.0°,  and  crosses  with  the  X  2B,  at  a  similar  angle.  Therefore  if  Of  is  produced  in  a  highly  bent 
geometry,  2B2  might  be  preferentially  populated.  The  large  transition  dipole  moment  and  the  vertical 

excitation  energy  for  2B2^2A,  also  suggest  that  2B2  may  be  electronically  excited  efficiently  to  2A,  at  the 

wave  length  of  523nm.  The  computed  energetics  of  2B2  and  2B,  and  the  corresponding  dissociation  limits 
may  explain  the  larger  maximum  KER  observed  in  the  second  experiment  of  Continetti  and  the  smaller 
02-0'  bond  energy  derived  from  the  experiment  of  Hiller,  if  we  assume  that  2B2  is  the  parent  state  in  both 
cases.  Furthermore,  meta-IRC  calculations  suggest  rather  different  final  state  distribution  of  the 

photofragments  from  2B,-^2A2  and  2B2^2A,  processes,  in  qualitative  agreement  with  the  experimental 

observations.  Although  the  vibrationally  excited  ground  state  Of  might  also  produce  rotational  hot, 
vibrational  cold  photofragments  through  the  angular  dependence  of  the  seam  between  the  two  diabatic 
excited  2A”states,  the  exact  effect  of  parent  vibrational  excitation  requires  future  dynamics  calculations.  At 
the  current  stage,  our  calculations  strongly  support  that  the  2B2  electronic  state  has  been  accessed  in  the 
second  experiment  of  Continetti  et  al. 

B-2.  Structure,  spectroscopy  and  reaction  of  (N20:H20)+  cation  cluster. 

J.  E.  Stevens,  M.  C.  Holthausen,  and  K.  Morokuma,  The  Potential  Energy  Surface  of  the  (N20  H20)+ 
Cluster  Cation,  J.  Chem.  Phys.,  Ill,  7766-7773  (1999). 

The  ground  state  potential  energy  surface  of  the  (N20-H20)+  cluster  cation  is  explored  with  ab 
initio  and  density  functional  theories.  B3LYP  optimizations  are  used  to  determine  the  structure  of  the 
products  of  the  dissociation  of  the  cluster  ions  as  well  as  possible  structures  for  the  clusters  themselves 
and  transition  states  that  connect  various  minima.  Energetics  for  all  optimized  structures  are  determined 
with  the  G2M(RCC,MP2)  method.  The  results  are  used  to  interpret  collision-induced  dissociation  (CID) 
experiments  which  study  the  cluster  ion,  and  which  find  that  the  cluster  dissociates  to  form  H20+  +  N20, 
N2OH+  +  OH,  and  N20+  +  H,0  products.  The  calculations  an  (N20-0H2)+  complex  as  well  as  a  similar 
(H20-N20)+  complex,  and  show  that  these  complexes  access  the  experimentally  observed  HzO+  +  N20 
products  and  N2OH+  +  OH  products  without  any  intervening  reverse  barrier.  The  stability  of  both  these 
complexes,  approximately  -20  kcal/mol  relative  to  the  H20+  +  N20  products,  agrees  well  with 
experimentally  determined  CID  thresholds  for  all  products.  Additional  calculations  of  the  ground  state 
potential  energy  surface  of  the  cluster  investigate  the  possibility  of  the  formation  of  other  products.  Some 
preliminary  studies  of  the  excited  states  of  the  cluster  cation  are  also  performed;  the  results  of  these 
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calculations  lend  insight  into  experimental  photodissociation  studies  of  the  cluster  ions.  Mechanisms  for 
the  formation  of  H20+  +  N20,  N2OH+  +  OH,  and  N20+  +  H20  products  following  photoexcitation  of  the 
cluster  ions  are  discussed;  the  H20+  +  N20  and  N2OH+  +  OH  products  must  be  formed  from  a  surface¬ 
hopping  from  an  excited  electronic  state  to  states  which  correlate  to  ground  state  products.  Similarly, 
N20+  +  H20  products  may  be  formed  from  collision  induced  dissociation  of  clusters  only  by  means  of  a 
surface-hopping  mechanism. 

B-3.  Electronic  structure  of  triply  charged  molecular  ions  N2+3. 

D.  Bandrauk,  D.  G.  Musaev  and  K.  Morokuma,  Electronic  States  of  the  Triply  Charged  Molecular  Ion 
N23+  and  Laser  Induced  Coulomb  Explosion,  Phys.  Rev.  A,  59,  4309-4315  (1999). 

Potential-energy  curves  for  the  low-lying  doublet  and  quartet  states  of  various  symmetries  (FIg,  II  u, 

Zg,  and  £u)  of  the  triply  charged  N23+  mol.  cation  are  studied  using  quantum  ab  initio  methods.  These 

states  correlate  asymptotically  to  the  N2+  +  N+  at.  fragments.  No  stable  (bound)  states  are  found  to  occur 
in  the  midst  of  many  Coulomb  repulsive  potentials.  The  potentials  are  shown  to  follow  a  pure  Coulomb 

repulsion  law  at  R>5  bohr.  The  most  relevant  configurations  are  enumerated  for  discussion  of  recent 
laser  Coulomb  explosion  and  imaging  experiments. 

B-4.  Structure  and  spectroscopy  of  Cl3‘  ion. 

A.  L.  Kaledin,  M.  C.  Heaven,  K.  Morokuma  and  D.  M.  Neumark,  Cl3'  electron  photodetachment 
spectrum:  measurement  and  assignment,  Chem.  Phys.  Lett.,  306,  48-52  (1999). 

An  electron  photodetachment  spectrum  of  Cl3~  is  reported  for  193  nm  (6.42eV)  excitation.  The 
spectrum  was  assigned  using  high-level  ab  initio  calculation  for  the  Cl3  radical  and  the  Cl3~  ion.  A  broad 
band  centered  around  1.25  eV  (electron  kinetic  energy)  has  been  assigned  to  the  Cl3“  (X  ’Z+g)  — >  Cl3  (X 
2nu)  and  Cl3  (1 2Z+g)  transitions.  Vertical  photodetachment  accesses  the  transition  state  region  for  the  Cl 
exchange  reaction:  Cl  +  Cl2  — >  Cl2  +  Cl.  A  narrow  doublet  band  at  0.55  eV  is  assigned  to  the  Cl3“  (X 
lX+g)  — >  Cl3  (1 2ng)  transition,  split  by  spin-orbit  interaction. 

B-5.  The  photoelectron  spectroscopy  of  the  IHBr  and  IDBr  anions. 

A.  Kaledin,  S.  Skokov,  J.  M.  Bowman  and  K.  Morokuma’  Theoretical  study  of  the  photoelectron 
spectroscopy  of  the  IHBr  and  IDBr  anions,  J.  Chem.  Phys.,  113,  9479-9487  (2000). 

The  coupled  cluster  method  with  a  large  basis  set  with  quasi-relativistic  effective  core  potentials 

on  the  halogens  was  utilized  to  investigate  the  ground  electronic  state,  X'Z+,  of  the  IHBr  anion.  A  semi- 


9 


global,  three-dimensional  potential  energy  surface  was  obtained  via  spline  interpolation  of  ab  initio  data. 
Variational  rovibrational  calculations  were  carried  out  with  this  potential  energy  surface,  and  assignment 
of  fundamentals  and  low  lying  overtones  as  well  as  combination  bands  is  presented  for  IH(D)Br  anions. 

Photodetachment  spectra  for  IH(D)Br  +  hv  -»  [IH(D)Br]  +  e  '  were  computed  using  three-dimensional 

time  independent  real  and  complex  L2  methods  and  a  recently  developed  real  L2  wave  packet  propagation 
technique;  all  methods  employed  a  previous  London-Eyring-Polanyi-Sato  empirical  potential  for  the 
neutral  system.  Well  resolved  fine  structure  is  found  in  photodetachment  spectra  of  IHBr.  This  structure 
is  assigned  to  bending  excitation  based  on  analysis  of  the  wavefunctions  and  the  neutral  potential  energy 
surface  near  the  Franck-Condon  region.  This  fine  structure  is  largely  eliminated  in  the  EDBr  spectrum. 
The  photodetachment  spectrum  for  the  first  excited  asymmetric  stretch  of  IHBr'  shows  a  new  feature  that 
is  associated  with  sampling  of  the  transition  state  region  of  the  neutral  surface.  Spin-orbit  corrected  multi¬ 
reference  configuration  interaction  calculations  indicate  that  the  first  excited  electronic  state  is  only  about 
0.06  eV  above  the  ground  electronic  state  of  the  neutral  IHBr  system  in  the  Franck-Condon  region.  The 
model  empirical  potential  used  in  the  calculations  of  photodetachment  spectra  is  shown  to  have  roughly 
the  average  behavior  of  these  two  ab  initio  potentials. 

C.  Potential  Energy  Surfaces  for  Photochemical  Reactions  and  Nonadiabatic  Processes. 
C-l.  Mechanism  of  photodissociation  of  isocyanic  acid,  HNCO. 

J.  E.  Stevens,  Q.  Cui,  and  K.  Morokuma,  An  ab  initio  Study  of  the  Dissociation  of  HNCO  in  the  S, 
Electronic  State,  J.  Chem.  Phys.,  108,  1452-1458  (1998). 

Regions  of  the  S,  potential  energy  surface  of  HNCO  relevant  to  N-H  and  C-N  bond 
photodissociation  have  been  investigated  with  ab  initio  calculations.  Geometries  of  minima  and  transition 
states  on  S,  as  well  as  those  of  the  product  photofragments  and  the  HNCO  ground  state  have  been 
optimized  with  the  CASSCF  method,  and  their  energies  calculated  with  MRSDCI  and  CASPT2  methods. 
Deep  planar  trans  and  cis  minima  exist  on  the  S,  surface,  and  are  connected  by  transition  states  for 

isomerization.  The  S0-»S,  electronic  transition  is  brighter  for  trans  configurations  than  for  cis,  and  the 

initial  excitation  and  dynamics  is  most  likely  to  proceed  through  trans  configurations.  The  N-H  fission  on 
S,  has  a  substantial  barrier;  it  occurs  more  easily  through  the  planar  cis  transition  state  which  is  about  20 
kcal/mol  above  the  dissociation  threshold  than  through  the  trans  transition  state.  The  C-N  fission  on  S[ 
can  take  place  both  via  the  planar  trans  and  the  planar  cis  transition  state  with  a  low  barrier  over  the 
dissociation  threshold;  the  reverse  barrier  is  estimated  to  be  a  few  kcal/mol. 

A.  L.  Kaledin,  Q.  Cui,  M.  C.  Heaven  and  K.  Morokuma,  Ab  initio  Theoretical  Studies  on 
Photodissociation  of  HNCO  upon  S,(‘A”)  <-  S0('A’)  Excitation:  The  Role  of  Internal  Conversion  and 
Intersystem  Crossing,  J.  Chem.  Phys.,  Ill,  5004-5016  (1999). 
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Photodissociation  of  isocyanic  acid,  HNCO,  was  studied  with  high-level  ab  initio  methods. 
Geometry  optimizations  of  stationary  points  and  surface  crossing  seams  were  performed  with  the 
complete  active  space  self-consistent  field  (CASSCF)  method,  and  the  energetics  were  re-evaluated  with 
single-point  second-order  multireference  perturbation  theory  (CASPT2).  The  three  product  channels  that 

participate  in  the  photodissociation  process  are:  [1]  HN(X3Z  )  +  CO  at  86.0  (calculated  79.6)  kcal/mol, 

[2]  H  +  NCO(X2n)  at  109.7  (108.7)  kcal/mol,  and  [3]  HN(a'A)  +  CO  at  122.2  (120.8)  kcal/mol.  The 

four  electronic  states,  S0,  S ,,  T,  and  T2,  that  interconnect  these  channels  were  studied  in  detail.  S, 
exhibits  dissociation  barriers  to  both,  channel  [2]  and  [3],  whose  respective  reverse  heights  are  11.3  and 
1.2  kcal/mol,  in  good  agreement  with  experiment  as  well  as  previous  theoretical  works.  The  two  triplets, 
T,  and  T2,  show  barriers  of  similar  heights  for  HN  bond  fission,  while  S0  has  no  barriers  to  either 
channel.  Various  key  isomerization  transition  states  as  well  numerous  minima  on  the  seam  of  surface 
crossings  (MSX’s)  were  also  found.  At  photoexcitation  energies  near  channel  [3]  threshold,  products  to 

channel  [3]  are  likely  to  be  formed  via:  S,-»[3]  (if  enough  energy  in  excitation)  and  S ,  — > S0 — > [3] . 
Channel  [2]  can  be  formed  via:  S,— >S0— »[2];  (HN-mode  quanta)+S,— »T,— »[2];  S,— »T2— >[2]; 
Sj — »T2 — >T[ — >[2),  and  channel  [1]  via  Sj— >S0— >Tt — >[1],  Sj — >Tj — >[1]  and  S, — >T2 — >T( — >[1],  At 

higher  photoexcitation  energies  S,— ■ >  [3]  is  expected  to  be  dominant  while  channel  [2]  with  higher 
activation  energy  is  expected  to  drop  rapidly.  Also  addressed  are  such  important  issues  as  the  impact  of 
vibrationally  excited  HN  mode  on  channel  [2]  yield,  and  the  band  origin  of  the  S  ,< — S0  excitation 
spectrum. 

C-2.  Mechanism  of  spin-forbidden  reaction  CH(2P)  +  N2  — »  HCN  +  N(4S) 

Q.  Cui  and  K.  Morokuma,  The  spin-forbidden  reaction  CH(2P)  +  N2  — »  HCN  +  N(4S)  revisited.  I.  — 

Ab  initio  study  on  the  potential  energy  surfaces,  Theo.  Chem.  Acc.,  102,  127-133  (1999). 

High  level  ab  initio  electronic  structure  theories  have  been  applied  to  investigate  the  detailed 

reaction  mechanism  of  the  spin-forbidden  reaction  of  CH(2II)  +  N2  — >  HCN  +  N(4S).  The  G2M(RCC) 

calculations  provide  accurate  energetics  for  the  intermediates  and  transition  states  involved  in  the  reaction, 
whereas  the  B3LYP/6-311G(d,p)  method  overestimates  the  stability  of  some  intermediates  as  much  as 
~10kcal/mol.  A  few  new  structures  have  been  found  on  both  the  doublet  and  quartet  electronic  states, 
which  are  mainly  involved  in  the  dative  pathways.  However,  due  to  the  higher  energies  of  these 
structures,  the  dominant  mechanism  remains  to  be  the  one  involving  the  C2v  intersystem  crossing  step. 
The  C2v  MSX  structures  and  the  spin-orbit  coupling  between  the  doublet  and  quartet  electronic  states  are 
rather  close  to  those  found  in  previous  studies.  Vibrational  frequencies  orthogonal  to  the  norm  of  the 
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seam  have  been  calculated  at  the  MSX  from  the  first  principle,  which  have  been  applied  in  a  separate 
publication  to  calculate  the  rate  of  the  titled  reaction  with  a  newly  proposed  non-adiabatic  transition  state 
theory  for  spin-forbidden  reactions. 

Q.  Cui,  K.  Morokuma,  J.  M.  Bowman  and  S.  J.  Klippenstein,  The  spin-forbidden  reaction  CH(2P)  +  N2 
— >  HCN  +  N(4S)  revisited.  II.  —  Non-adiabatic  transition  state  theory  and  application.  J.  Chem.  Phys., 
110,  9469-9482  (1999). 

Transition  state  theory  is  extended  straightforwardly  to  treat  non-adiabatic  processes  and  applied 
to  study  the  rate  constant  of  the  spin-forbidden  reaction  CH(2II)  +  N2  —>  HCN  +N(4S).  A  one¬ 
dimensional  model  was  set  up  to  calculate  the  intersystem  crossing  probability  with  the  distorted  wave 
approximation  and  using  an  ab  initio  value  of  the  spin-orbit  coupling.  The  effect  of  orthogonal  degrees  of 
freedom  was  then  considered  by  energy  convolution  with  the  vibrational  frequencies,  obtained  from  ab 
initio  calculations,  orthogonal  to  the  crossing  seam  at  the  minimum  of  the  seam  of  crossing  (MSX),  also 
obtained  from  ab  initio  calculations.  An  expression  for  the  cumulative  reaction  probability,  N(E),  of  the 
reaction  was  obtained  by  a  straightforward  extension  of  the  unified  statistical  theory,  where  the  MSX  was 
treated  as  a  transition  state.  The  calculated  thermal  rate  constant,  k(T),  seems  to  be  too  low  by  two  orders 
of  magnitude  compared  to  experimental  measurements  and  an  empirical  transition  state  study  where 
empirical  vibrational  frequencies  at  the  MSX  are  lower  by  a  factor  of  two  than  those  derived  here.  The 
disagreement  strongly  suggests  that  the  current  treatment  of  the  multi-dimensional  dynamics  needs  to  be 
improved.  In  particular,  it  may  be  a  poor  assumption  that  spin-forbidden  transition  takes  place  with 
uniform  probability  on  the  seam  in  the  case  we  are  considering. 

C-3.  Ab  initio  direct-trajectory  study  of  the  photodissociation  of  CIOOC1. 

L.  Kaledin,  and  K.Morokuma,  An  ab  initio  direct- trajectory  study  of  the  photodissociation  of  ClOOCl,  J. 
Chem.  Phys.,  113,  5750-5762  (2000). 

The  photodissociation  of  chlorine  peroxide,  ClOOCl,  is  studied  with  classical  trajectories  where 
the  energy  and  gradient  are  computed  on  the  fly  by  means  of  the  state-averaged  (sa)  CASSCF  with  the 
DZP(+)  basis  set.  We  show  that  six  electronically  excited  states  are  involved  in  the  process  of 
decomposition,  which  proceeds  via  several  competing  pathways  and  at  least  three  electronically  unique 
fragment  channels.  The  problem  is  treated  in  4D  (C2  constraint)  and  5D  (planar  constraint)  frameworks  in 
order  to  model  the  mechanisms  of  synchronous  and  asynchronous  or  stepwise  dissociation,  respectively. 
A  single  trajectory  with  the  initial  conditions  of  a  non-vibrating  non-rotating  molecule  is  propagated  on 
each  excited  state  surface  for  an  average  time  of  10  fs  for  the  purposes  of  determining  the  early  stages  of 
bond  breaking.  We  show  that  even  in  such  a  short  propagation  time  the  pathway  competition  can  be  more 
or  less  unambiguously  understood.  The  results  indicate  that  in  the  regime  of  a  308  nm  photolysis,  the 
major  dissociation  fragments  are  Cl  atoms  and  02  molecules,  both  in  the  ground  state.  The  higher  energy 
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regime  of  a  248  nm  photoexcitation  yields  additional  fragments,  e.g.,  C10(X2n),  0(3P)  and  C100(X2A”, 

12A’).  We  have  achieved  an  overall  qualitative  agreement  with  experiment  that  more  than  70%  of  the 
available  energy  is  transferred  into  the  translational  energy  of  the  products  for  the  case  of  the 
synchronous  concerted  dissociation.  In  all  the  cases,  the  rotational  excitation  of  produced  molecular 
oxygen  is  very  high,  while  its  vibration  is  in  v=0.  Implications  of  the  results  on  the  stratospheric  ozone 
depletion  cycle  are  also  presented. 

C-4.  An  initio  direct-trajectory  study  of  the  photodissociation  of  ketene. 

A.  L.  Kaledin,  J.  Seong  and  K.  Morokuma,  Predominance  of  Non-equilibrium  Dynamics  in  the 
Photodissociation  of  Ketene  in  the  Triplet  State,  J.  Phys.  Chem.  A,  105,  2731-2737  (2001). 

The  photodissociation  of  ketene  is  studied  using  direct  surface  hopping  classical  trajectories 
where  the  energy  and  gradient  are  computed  on  the  fly  by  means  of  state-averaged  complete  active  space 
self-consistent  field  with  a  double  zeta  polarized  basis  set.  Three  low-lying  electronic  states,  singlets  S0, 
S,  and  triplet  T,,  are  involved  in  the  process  of  photodissociation  of  triplet  state  ketene.  We  propagated  a 
trajectory,  starting  at  the  Franck-Condon  geometry  on  S,,  and  branched  it  out  into  many  child  trajectories 
every  time  the  propagating  potential  energy  surface  (PES)  crossed  with  another  PES.  The  major 

photodissociation  pathway  to  the  triplet  products  was  found  to  be  S,  — >  S0  — >  Tj  — »  CH2(X  3B1)  +  CO(X 

’E+).  It  has  been  found  that  1)  the  S0-T,  non-adiabatic  transition  creates  the  T,  species  non-statistically  at 

restricted  regions  of  phase  space,  and  2)  a  large  fraction  of  the  T,  species  thus  created  dissociates  almost 
immediately,  leaving  no  time  for  equilibration  of  internal  degrees  of  freedom.  Whether  a  specific  T, 
trajectory  dissociates  fast  or  not  is  determined  by  the  amount  of  C-C  stretch  vibration  at  the  S0-T,  branch 
point.  In  essence,  the  above  observations  suggest  strongly  that  the  T,  photodissociation  process  is  highly 
non-statistical,  thus  making  equilibrium-based  statistical  theories  inapplicable  for  computing  the 
dissociation  rate. 

C-5.  The  symmetry  breaking  and  its  effect  on  the  potential  surface  of  N03. 

W.  Eisfeld  and  K.  Morokuma,  A  detailed  study  on  the  symmetry  breaking  and  its  effect  on  the  potential 
surface  of  NOv  J.  Chem.  Phys.,  113,  5587-5597  (2000). 

The  tenacious  symmetry  breaking  of  the  electronic  wave  function  of  the  nitrate  radical  (N03)  and 
its  effect  on  the  ground-state  potential  energy  surface  is  investigated  in  detail.  The  symmetry  breaking  of 
Hartree-Fock  wave  functions  results  from  a  dominance  of  the  orbital  localization  effect  over  the 
resonance  effect  and  leads  to  three  different  solutions,  one  symmetrical  and  two  distorted  ones,  for  the 
same  electronic  state.  The  respective  equilibrium  geometries  of  these  solutions  are  points  on  different 
potential  surfaces,  making  their  comparison  meaningless.  The  resonance  effect  is  promoted  by  dynamic 
as  well  as  static  electron  correlation.  However,  the  dynamic  correlation  methods  [e.g.,  many-body 
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perturbation  theory  (MBPT)  and  coupled-cluster  single  double  (CCSD)]  cannot  overcome  the  symmetry 
breaking  of  the  reference  function  and  the  problem  of  multiple  solutions  persists.  The  symmetry  breaking 
can  be  avoided  by  the  complete  active  space  self-consistent  field  (CASSCF)  approach  that  yields  unique, 
single-valued  surfaces  for  all  electronic  states.  However,  a  sufficiently  large  and  appropriately  selected 
active  space  has  to  be  used  to  avoid  unphysical  distortion  of  the  wave  function.  Still  the  orbital 
localization  effect  leads  to  equilibrium  geometries  of  C2v  symmetry  which  strongly  depend  on  the  state¬ 
averaging  of  the  CASSCF  wave  function.  Multireference  single  double  configuration  interaction  (MR- 
SDCI)  wave  functions  are  also  free  of  symmetry  breaking,  if  the  reference  orbitals  are  and  if  the 
configuration  space  is  invariant  under  the  symmetry  operations.  MRCI  geometry  optimizations  only 
result  in  D3h  symmetric  structures  with  bond  lengths  and  harmonic  frequencies  in  close  agreement  with 
experimental  data. 

W.  Eisfeld  and  K.  Morokuma,  Ab  initio  investigation  of  the  vertical  and  adiabatic  excitation  spectra  of 
N03,  J.  Chem.  Phys.,  114,  9430-9440  (2001). 

The  ground  and  excited  electronic  states  up  to  approximately  100  000  cm'1  of  the  nitrate  free 
radical  (N03)  are  investigated  by  complete  active  space  self-consistent-field  (CASSCF)  and 
mulitreference-singles  doubles  configuration  interaction  calculations.  Extended  basis  sets,  containing 
polarization  and  diffuse  functions,  and  an  active  space  consisting  of  13  orbitals  and  17  electrons  are  used. 
A  total  of  28  electronic  states  is  obtained  within  the  D3h  framework  which  split  into  44  states  in  C2v  All 
calculated  states  are  valence  states  and  their  character  is  discussed  in  detail.  Oscillator  strength  and 
radiative  lifetimes  are  determined  from  the  CASSCF  wave  functions  and  give  evidence  for  a  strong 
transition  that  was  not  yet  known  experimentally.  For  the  experimentally  observed  2E  and  2E  states 
equilibrium  geometries,  harmonic  frequencies,  and  adiabatic  excitation  energies  are  calculated  in  excellent 
agreement  with  experimental  data.  Important  new  insight  is  gained  about  the  role  of  the  low-lying 
electronic  states  in  the  photodissociation  and  the  mechanism  of  this  process  is  discussed. 

C-6.  Spin-allowed  and  Spin-forbidden  Pathways  of  the  Gas  Phase  Reactions  of  0(3P) 
+  C2H5I, 

J.  E.  Stevens,  Q.Cui  and  K.  Morokuma,  An  ab  initio  Investigation  of  Spin-allowed  and  Spin-forbidden 
Pathways  of  the  Gas  Phase  Reactions  of  0(3P)  +  C2H5I,  J.  Chem.  Phys.,  108,  1544-1551  (1998). 

The  singlet  and  triplet  potential  energy  surfaces  involved  in  the  gas  phase  reactive  collisions  of 
0(3P)  and  C2H5I  have  been  studied  with  ab  initio  electronic  structure  computations.  The  collisions 
produce  both  spin-forbidden  HOI  +  C2H4  and  spin-allowed  OI  +  C2H5  products.  The  calculations 
indicate  that  HOI  is  formed  via  a  triplet  complex  and  through  a  triplet/singlet  intersystem  crossing, 
followed  by  passage  through  a  singlet  intermediate  and  transition  state  for  the  intramolecular  abstraction 

of  P-hydrogen.  All  the  relevant  structures  for  this  pathway  are  lower  in  energy  than  the  reactants,  and  this 
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pathway  is  accessible  even  at  low  impact  energies.  The  calculations  also  indicate  that  01  may  be  formed 
by  two  channels.  One  is  the  same  to  the  above  singlet  pathway  up  to  the  singlet  intermediate,  which  now 
dissociates  endothermically  without  barrier  to  give  the  products.  The  second  channel  is  the  direct 
dissociation  of  the  triplet  intermediate,  and  is  open  only  when  an  enough  excess  energy  to  surmount  a 
triplet  transition  state  is  provided.  The  product  energy  distribution  is  also  discussed  based  on  the 
structures  of  transition  states. 

D.  Potential  Energy  Surfaces  of  Open-Shell  Molecules. 

D-l.  Electronic  structure  of  C2H. 

Q.  Cui  and  K.  Morokuma,  Ab  initio  Studies  on  the  Electronic  Excited  States  of  C2H,  J.  Chem.  Phys., 
108,  626-636  (1998). 

High  level  ab  initio  calculations  have  been  carried  out  to  study  seven  electronic  states  of  C2H.  The 
calculated  equilibrium  structure,  energetics  and  vibrational  frequencies  for  the  3  2A’  state  at  the 
CASPT2/PVTZ  level  are  in  good  agreement  with  those  obtained  experimentally  by  Hsu  et  al.  The 
transition  dipole  moments  from  the  ground  state  C2H  to  electronic  excited  states  depend  sensitively  on  the 
H-C-C  bending  angle  and  often  peak  at  non-linear  configurations.  Based  on  this  and  the  dissociation 
behavior  of  the  excited  states,  we  predict  A  'IIU  and  c  3Z+U  C2  fragments  to  be  rich  in  population,  the 
former  of  which  is  experimentally  detected  recently  by  Jackson  et  al.  The  ground  X  'Z+g  and  the  a  3I1U 

state  C2  are  expected  to  be  formed  via  the  non-adiabatic  process  3  2A’  — >  2  2A’  or  4  2A’  — >  3  2A’  — >  2 

2A\  which  is  in  accord  with  the  experimentally  observed  lifetime  pattern  by  Hsu  et  al.  No  reverse  barrier 
for  CC-H  dissociation  was  found  on  the  X  and  A  electronic  states  of  C2H  in  the  linear  configuration.  The 
2  2 A'  state,  however,  develops  a  distinct  “barrier”  (not  a  true  saddle  point)  along  dissociation  coordinate 
when  the  H-C-C  is  significantly  bent,  due  to  the  interaction  with  upper  electronic  state.  Since  the  2  2A’ 
state  energetically  prefers  a  linear  dissociation,  we  suspect  that  the  upper  bound  of  the  CC-H  bond  energy 
measured  by  Hsu  et  al.  is  not  severally  affected  by  this  “barrier”. 

D-2.  Structure,  electronic  structure  and  spectroscopy  of  Cl3. 

A.  L.  Kaledin,  M.  C.  Heaven,  W.  G.  Lawrence,  Q.  Cui,  J.  E.  Stevens  and  K.  Morokuma,  Ab  initio 
Molecuclar  Orbital  Study  of  the  Trichlorine  Radical,  Cl3,  J.  Chem.  Phys.,  108,  2771-2783  (1998). 

We  report  a  rigorous  ab  initio  study  of  the  ground  and  low-lying  excited-state  potential  energy 
surfaces  (PESs)  of  the  Cl3  radical  at  CASSCF,  CASPT2  and  MRSDCI  levels  of  theory.  The  ground  state 

has  two  C1..C12  van  der  Waals  complexes,  XL  and  X'B .  The  linear  asymmetric  minimum  (XL)  is  2n, 
with  a  Cl-Cl  distance  r  =  3.90  bohr,  and  a  Cl-M  (M:  the  Cl2  center-of-mass)  distance  R  =  8.70  bohr.  The 
bent  asymmetric  minimum  (X'B)  is  of  2A’  symmetry,  with  r  =  3.90  bohr,  R  =  6.85  bohr,  and  the  angle 

A  ~  _ 

between  r  and  R,  y=  68.4°.  Spin-orbit  Cl  predicts  that  the  global  minimum  is  linear  XL  (  fI3/2)  with  a 
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bond  dissociation  energy  of  De(Cl2(X)-Cl)  of  280  cm1.  Low-lying  doublet  excited  states  have  only  one 
strongly  bound  structure,  a  linear  symmetric  AL  (1  2llg)  state  with  a  bond  distance  of  4.67  bohr.  This 

state  is  bound  by  approximately  4,300  cm'1  with  respect  to  the  Cl2(3riu)  +  Cl  asymptote,  and  its  minimum 

lies  about  8,700  cm'1  above  the  XL  van  der  Waals  minimum.  Transition  dipole  moment  calculations 

show  that  the  A  —  X  transition  is  fully  allowed.  Two  bound  quartet  minima  were  located.  The  most 
deeply  bound  was  QD3h  (14A,')  with  a  D3h  equilibrium  geometry  (r  =  5.00  bohr)  about  11,300  cm'1 

above  XL.  The  other  state,  QC2v  (14A2)  had  a  C2v  equilibrium  geometry  (r,  =  4.83  bohr  and  6  = 

101.7°)  and  an  energy  of  about  13,500  cm'1  relative  to  XL.  Although  C13(X)  is  shown  to  be  unstable, 
the  present  results  support  the  notion  that  Cl3  participates  in  Cl  atom  recombination  processes.  However, 
the  energies  and  transition  moments  of  the  low-lying  excited  states  are  not  consistent  with  electronic 
spectra  that  have  been  tentatively  assigned  to  Cl3 

E.  Dynamics  and  Kinetics  in  the  Ground  Electronic  State. 

E-l.  H  and  H2  elimination  pathways  from  methane,  ethane  and  propane. 

S.  Irle  and  K.  Morokuma,  A  molecular  orbital  study  on  H  and  H2  elimination  pathways  from  methane, 
ethane  and  propane,  J.  Chem.  Phys.,  113,  6139-6148  (2000). 

Decomposition  pathways  for  propane,  as  well  as  methane  and  ethane  for  comparison,  in  its 
ground  electronic  state  were  studied  using  density  functional  and  high  accuracy  ab  initio  calculations.  The 
reaction  pathways  were  characterized  by  locating  the  transition  states  and  following  the  intrinsic  reaction 
coordinate.  Atomic  hydrogen  as  well  as  molecular  hydrogen  elimination  pathways  were  investigated, 
including  three  deuterated  propane  species  for  comparison  with  experiment.  The  methyl  and  ethyl  groups 
in  propane  are  found  to  stabilize  transition  states  and  radical/biradical  intermediates  along  the  reaction 
pathways.  For  propane,  2,2-elimination  of  an  hydrogen  molecule  is  found  to  be  the  dominat  molecular 
elimination  pathway,  in  agreement  with  recent  photochemical  experiments.  We  find  1,1 -elimination  to  be 
the  next  important  molecular  elimination  pathway,  followed  by  1,2-elimination,  which  disagrees  with  the 
experimental  result  favoring  1,2-  over  1,1 -elimination. 

E-2.  Dipole  moments  of  highly  vibrationally  excited  HCN. 

J.  M.  Bowman,  S.  Irle,  K.  Morokuma,  and  A.  Wodtke,  Dipole  moments  of  highly  vibrationally  excited 
HCN:  theoretical  prediction  of  an  experimental  diagnostic  for  delocalized  states.  J.  Chem.  Phys.,  114, 
7923-7934  (2001). 

Vibrational  state  specific  dipole  moments  are  diagnostic  of  the  degree  of  localization  of  vibrational 
states  in  highly  vibrationally  excited  HCN.  Using  a  newly  calculated  global  ab  initio  dipole  moment 
function  and  previously  calculated  highly  accurate  vibrational  wave  functions,  we  show  that  delocalized 
(i.e.  isomerizing)  vibrational  states  of  HCN  possess  markedly  lower  dipole  moments  than  localized  HCN 
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or  HNC  states.  We  also  show  that  the  vibrational  quantum  number  dependence  of  the  dipole  moment  can 
be  used  to  distinguish  delocalized  states  from  localized  Franck-Condon-dark  states  that  are  made 
observable  by  perturbations  with  localized  Franck-Condon-bright  states.  Furthermore,  using  classical 
trajectory  analysis  we  introduce  and  describe  a  new  experimental  approach  to  obtain  these  data,  which 
relies  on  combining  optical  pumping  and  state  specific  molecular  transport  with  hexapoles.  With  this 
method  it  is  possible  to  determine  state  specific  dipole  moments  with  high  accuracy  and  precision. 

F.  Potential  Energy  Surfaces  for  the  Dynamics  of  Chemical  Laser  Systems. 

In  collaboration  with  Prof.  Mike  Heaven  and  the  Chemical  Laser  group  at  the  Air  Force  Research 
Laboratory,  Kirtland  Air  Force  Base,  we  have  undertaken  computational  studies  on  the  potential  energy 
surfaces  and  dynamics  of  chemical  reactions  and  energy  transfer  processes  which  control  the  efficiencies 
of  chemical  laser  systems.  In  particular,  we  have  worked  on  the  I  +  02  and  the  I +NC1  systems. 

F-l.  A  surface-hopping  trajectory  study  for  theoretical  prediction  of  the  rate  constant 
for  the  I(2P3/2)  +  C^fa’Ag)  <=>  I(2P1/2)  +  02(X3X'g)  energy  transfer  process. 

A.  L.  Kaledin,  M.  C.  Heaven,  and  K.  Morokuma,  Ab  initio  Potential  Energy  Surfaces  for  the  I(2P3/2)  + 
02(a'Ag)  o  I(2P|/2)  +  02(X3ITg)  Energy  Transfer  Process,  Chem.  Phys.  Lett.,  289,  110-117  (1998). 

Ab  initio  electronic  structure  calculations  have  been  used  to  investigate  the  electronic  energy 
transfer  process 

I(2P  3/2)+02(a‘ Ag)<=>I(2P  1/2)+02(X3E"g) 

Potential  energy  surfaces  for  all  states  associated  with  the  reactants  and  products  were  obtained  using 
CASSCF  and  CASPT2  methods,  including  the  effective  one-electron  spin-orbit  Hamiltonian.  Surfaces 
correlating  with  the  above  reactants  and  products  were  all  found  to  be  non-bonding.  Shallow  van  der 
Waals  minima  were  predicted  at  long  range.  Surface  crossings  were  found  at  energies  below  the 

I(2PV2)+02(alAg)  asymptote.  It  is  probable  that  these  crossings  are  responsible  for  the  efficient  transfer  of 

electronic  energy  in  this  system. 

L.  Kaledin,  M.  C.  Heaven,  K.  Morokuma,  Theoretical  prediction  of  the  rate  constant  for  I  +  02(a‘Ag) 

electronic  energy  transfer:  a  surface-hopping  trajectory  study,  J.  Chem.  Phys.,  114,  215-224  (2001). 

The  temperature  dependence  of  the  rate  constant  for  the  electronic  energy  transfer  process  I(2P3/2) 

+  02(a'Ag)  — >  I(2P1/2)  +  02(X3Eg  )  has  been  studied  theoretically.  Seven  ab  initio  diabatic  potential  energy 

surfaces,  four  for  the  entrance  channel  and  three  for  the  exit  channel,  and  the  coupling  elements  between 
them  were  adopted.  Energy  transfer  dynamics  was  simulated  with  the  semi-classical  surface-hopping 
trajectory  calculation,  using  Tully’s  “fewest  switches”  model  for  electronic  transition.  Approximately 
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5xl05  trajectories  were  statistically  averaged  over  a  range  of  impact  parameters  and  collision  energies  to 
calculate  thermal  rate  constants  for  the  temperature  range  10-300  K.  It  was  found  that  collisions  resulting 
in  energy  transfer  were  dominated  by  single  hop  trajectories.  The  calculated  energy  transfer  rate  constant 
was  found  to  decrease  smoothly  with  increasing  temperature  over  the  range  100-300  K.  The  predicted 
value  was  in  excellent  agreement  with  the  experimental  result  for  150  K,  but  the  calculations 
underestimate  room  temperature  data  by  a  factor  of  1.6.  The  rate  constant  increases  with  decreasing 
energy  because  (i)  long-range  attractive  forces  draw  slow  moving  collision  partners  together  and,  (ii) 
longer  lifetime  of  slow  collisions  increases  the  probability  of  surface-hopping.  It  is  also  found  that  there 
is  a  competition  between  rotational  relaxation  of  02(a)  and  electronic  energy  transfer 
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